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1. INTRODUCTION 

Superconductivity is a macroscopic manifestation of quantum mechanics. This is notice- 
able in Josephson junctions and in the Little-Parks effect |p, where the behavior of a ring is 
■ not determined by the local magnetic field, but rather by the magnetic flux which it encloses. 
0> ■ In the original Little-Parks experiment, the superconducting sample was a long cylindrical 
bJO' shell, intended to have uniform thickness, embedded in a magnetic field parallel to the axis. 

. The sample 'does not like' to enclose a non-integer number of magnetic flux quanta. As a con- 
00 ' 

sequence, when the flux deviates from an integer number, the transition to superconductivity 
occurs at a lower temperature. 
^ i Another interesting feature of superconductivity is its globallity. More precisely, a uniform 

^ I material at uniform temperature will be superconducting at some places and normal at others 
Q^ ' if this nonuniform pattern minimizes the overall thermodynamic potential. When parts of the 
Q ■ sample are normal, while the rest is superconducting, the material is said to be in a mixed 

o ■ 
u 
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^ ■ ductors. In this case normal and superconducting domains of finite volume coexist. Another 
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state. A known case of mixed state is the "intermediate state" M, 11 1 for type I supercon 



known case is that of "vortices" 0, |rT|] for type II superconductors. For each vortex, the 



material is normal (i.e. the order parameter vanishes) along a line. Using physical arguments 
and numerical minimizations, we recently predicted a situation in which quantization and 
globallity give rise to a qualitatively different kind of mixed state. In this case the sample is 
doubly connected, i.e. homotopic to an annulus, but it is energetically advantageous to break 
the superconductivity at a suitable layer, so that the superconducting part is singly connected. 
When this is the situation, we shall say that the sample is in the "singly-connected phase" . 

In the present paper we analyze the different kinds of functions that minimize the ther- 
modynamic potential in the Little-Parks experiment. They depend upon the temperature 
and the magnetic flux, with the deviation from uniform thickness as an additional parameter. 



The deviation from uniformity will be usually taken as a perturbation. A complementary 
analysis is provided by the special case of shells whose thickness is piecewise constant. We 
shall see that the "phase diagram" (the boundaries between qualitatively different minimizers 
in the temperature-flux plane) contains several bifurcation lines and critical points. Even an 
infinitesimal deviation from uniform thickness turns out to produce a qualitative change in 
the phase diagram. We provide an analytic treatment of the minimization problem near the 
bifurcation lines and the critical points. 

In the next section we formulate our energy functional, define nondimensional quantities 
and derive the Euler-Lagrange equations. In Section 3 we consider the case of uniform rings 
and gain some insight into the phase transition picture. Bifurcation lines from the normal 
state are studied in Section 4. We show that there actually exist two types of bifurcations: 
one in which the sample becomes superconducting everywhere, and another which involves a 
mixed (singly-connected) state. The mixed state itself may bifurcate into doubly connected 
superconducting states. The location of such bifurcations in the parameters plane (magnetic 
flux and temperature) is found in Section 5. 

We show that there exist two critical points where bifurcation lines intersect. Section 6 
is devoted to a detailed analysis of the various phases near these points. In Section 7 we 
summarize our results and discuss them. In particular we propose several experimental setups 
to verify our theoretical predictions. 

2. FORMULATION 

We deal with a sample of superconducting material with the topology of a ring. The traditional 
treatment of the Little- Parks experiment may be found in Here we develop a formalism 



which is suited for the case of nonuniform thickness. The free energy of the sample can be 
expressed by means of the Ginzburg-Landau theory @, |@], HI, 0- If we write the order 
parameter in the form ip = ■ue*'^, then the free energy density of the superconducting phase, 
relative to the normal phase, is 

g = -au" + ]^f3u^ + (H - //ext)V(8^) + l^'iVtiY + u\hVcP - e* A/cf]/{2m*). (1) 



Here e*/2 and m*/2 are the electron charge and mass, H = V x A is the magnetic field, i^ext 
is the magnetic field which would be present in the absence of superconductivity, h is Planck's 
constant divided by 2tt, c is the speed of light, and are two material parameters. We 
denote by Tc the critical temperature in the sample without magnetic field and are interested 
in the region where the temperature T is lower than and close to T^. The parameter a is 
proportional to 1 — T/Tc. 

To fix ideas, we shall consider a sample in the form of a cylindrical shell and use cylindrical 
coordinates {z,r,6), with the axis of the coordinates at the axis of the shell. We denote by 
R the radius of the cylinder (measured up to the middle of the shell) and denote by D{6) 
the area obtained by cutting the shell with a plane of given 6. D{9) is the product of the 
"height" (parallel to the axis) and the "thickness" (along the radial direction). The height will 
be taken as uniform. In the original Little-Parks experiment, the height was much larger than 
the thickness; in recent experiments, the height is smaller than the thickness. There are two 
lengthscales in superconductivity: the magnetic field varies significantly over a lengthscale 
called the penetration depths given by \l = {c/ e*){m* (3 / AixaY^"^; the lengthscale over which 
the order parameter varies significantly is the coherence length ^ = jlrrt'a)^!'^ . Since 
a — * for T T^, both lengthscales are large near Tg. We shall see that the interesting 
temperatures are those where ^ is comparable with R. We shall take thickness and height 
much smaller than R, so that ip becomes a function of 6 only. (Also in the original Little-Parks 
experiment ip is independent of z, due to translational symmetry.) 

Since H — H^xt is proportional to the current along the ring, its contribution to the free 
energy is proportional to the square of D{6), whereas the contribution of the other terms is 
linear in D{6). Therefore, the term with if — ifgxt becomes negligible at the interesting temper- 
atures, and will be dropped for simplicity. The contribution of this term will be evaluated in 
Appendix A and, in particular, we shall have a criterion for the validity of this approximation. 

In order to write the free energy in a more convenient form, we observe that the supercur- 
rent I along the ring is given by 

I = {e*/m*)u\hV(j) - e*A/c)D{e), (2) 
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and it is independent of 9. Dividing (0) by u^D, and integrating along the ring, we obtain an 
alternative representation for the supercurrent, 

/ = 2Txkhae*/{m*(3K). (3) 

Here 

fc = ^[0(277) -</)(0)]-$/$o, (4) 
$ is the magnetix flux, $o = '^^T^Tic/e* is the flux quantum, 

A = R y-^D-^de, (5) 

JO 

and we have scaled the order parameter y = {(S/aY^'^u. 

We can now use the equality between both expressions for I to get rid of the term with 
V0 in (|l]) and write the energy functional in the form 

G=f ^= ri-y' + \y' + dnyyMe)de + (6) 

Jsample a^/t JO 2 K HA 

Note that in obtaining (P) we have not actually required any particular shape of the 
ring. We shall keep this notation with circular appearence for aesthetic reasons, but all our 
treatment remains valid for an arbitrary loop if we define R as its length divided by 2tt and 6 
as the arc length divided by R. Also for aesthetic reasons, we define A = {R/^Y, multiply G 
by the constant A/ J^^ D{6)d9, take R as the unit of length, and arrive at our model for the 
free energy 

G[i/,0] = ^ (-Xy' + -Xy^ + (y Y)Dde + i27rk)'A-'j / Dd9. (7) 

A rigorous derivation of as a limit of the Ginzburg-Landau functional in thin rings will 
be published elsewhere. 

While the dependence of G on the amplitude y is quite complex, G depends on the phase 
only through k"^. For a given flux $, the phase will adjust itself to minimize fc^ in (^. 

In order to minimize G, the order parameter has to choose between two distinct strategies, 
each of which gives rise to a constraint. 
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• Doubly-connected phase (DC): The entire sample is superconducting {y > for all 6) 
and the topology of the superconducting region is the same as that of the sample. The 
constraint is that [4>{27t) — 0(O)]/(27r) must be integer and this condition restricts the 
minimum value which k"^ can assume for a given $. For any $, k"^ can be in the interval 
[0, |], and therefore it will be in this interval when G is minimized. The interval k"^ <\ 
will be dubbed the relevant region. However, there is much to be learned by allowing 
values A;^ > i in our model functional G, and we shall study this situation, too. The last 
term in the numerator of (|^) may be regarded as the energy price that superconductivity 
has to pay for staying doubly-connected and enclosing a noninteger flux. 

• Singly- connected phase (SC): In this case y vanishes for some 9, so that the supercon- 
ducting region is singly-connected and does not enclose the magnetic flux. 0(27r) and 
0(0) are now unrelated, so that their difference can adjust itself to yield = in (H) and 
the last term in the numerator of (0) vanishes. We can also regard this vanishing as due 
to the divergence of A (and then the value of k is actually irrelevant). The constraint is 
now that y vanishes at some 9, and the energy price is in the first and third term in the 
integrand in (0). We shall often consider the case in which D is an increasing function 
of 9 for < < TT and is symmetric about 9 = Q and 9 = n. We shall denote such D^s 
as functions of type S. In this case, the place with the lowest energy price for y to vanish 
will be at 6* = 0. 

Our goal is now to find the function y that minimizes G for given values of the parameters 
A, which is controlled by the temperature, and k, which is governed by the magnetic flux. 
In particular, we want to find the values of A and k for which y{0) vanishes, and how y{0) 
approaches as these parameters vary. Note that the minimum value of k"^ which is consistent 
with (^) is a periodic function of the flux. This implies that all physical properties will be 
periodic functions of the flux. This is a fair idealisation of the experimentally known results. 

Minimizing (|^ over all periodic functions, we obtain the Euler-Lagrange (EL) equation 

{Dy')' + XD{y-y')-{^-^r^y-' = 0. (8) 
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The task of solving Eq. (P) is complicated by the presence of the nonlocal quantity A. This 
reflects the nonlocal nature of the effect we are discussing. For the SC phase the EL equation 
reduces to 

{Dy')' + \D{y-y^) = Q, ^(0) = y(27r) = 0. (9) 

It is sometimes more convenient to work with w = y"^. We therefore write the functional 
G and the associated EL equation in terms of w 

G[w] = ( r{-Xw + Ixw^ + ^)Dde + I BdQ, (10) 

yjo 2 4w J\ j Jo 

Dww" + D'ww - ^{w'Y + 2\Dw\l - w) - = 0, (11) 

and now A = {wD)^^d9. If w vanishes at some point, then, as before, the terms involving 
A are not present. Adding to the numerator of (p!OD half the integral of (pJ]) divided by w, we 
obtain that, if the EL equation is fulfilled, then the free energy is 

Gel = --A ^ Dw^dO/ f^" Dd9 . (12) 
2 JO Jo 

This shows that any nontrivial solution of the EL equation has lower free energy than the 
normal state. 

Since G is invariant under multiplication of D{6) by a constant, the average of D can be 
set as 1 without loss of generality. Except for Subsection 4.3, we shall adopt this convention 
in the following. 

3. UNIFORM D AND OUTLINE OF GENERAL D 

The case i5=constant is the traditional Little-Parks problem. In this case it is not difficult 
to guess the solutions of (P) and (||). For the doubly-connected phase we try the ansatz 
2/= const ant, which gives 

w = y^ = 1- k'^/X, (13) 

for which G becomes 

Gconst = ~{\-ky/X. (14) 
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For the SC phase 



where m obeys 



this gives the free energy 



.sn(J-^^|m), (15) 



1 + m \ 1 + m 



K{m) = nJ-^; (16) 
V 1 + m 



2A , 2 + m mm) 
6{l + m) 1 + m K[m) 

In the equations above, sn is a Jacobian elliptic function and K and E are complete elliptic 
integrals [Q]. 

For the constant phase, we note that ?/ = for k"^ = X and becomes meaningless for k"^ > A. 
This defines a bifurcation hne from the normal to the constant phase, denoted by Fi in Fig. 1. 
Likewise, the SC solution vanishes for A = | and is meaningless for A < |. This defines the 
bifurcation line Fn from the normal to the SC phase. For k = 0, Gconst < Csc, but Gconst 
increases with k^. Gconst = Gsc for 

k' = \ + {l- - i) - - ^)Vv^ + 5(A - i)V54^/^ + ... (18) 

Except for the point P = (/c^ = A = where Gconst = Csc = 0, equation ([T8|) describes a 
line in the region k"^ > j. Therefore, the global minimizer in the relevant region is the constant 
phase. 

Let us check now whether there are bifurcations from the constant phase. For this purpose, 
let us anticipate the nonuniform case and write 

D = l + eDi{e) (19) 

where e is a small positive number and the deviation from uniformity is defined such that 
Jq^ Di{6)d6 = 0. Expanding the order parameter in the form 

w = l-k'^/X + ewi + ..., (20) 

with A > fc^, Eq. (|lTD gives to first order in e a linear nonlocal differential equation for wi. 
The nonlocal term vanishes by periodicity, and we obtain 

w'l + (6fc2 - 2X)wi = -4^(1 - k^/X)Di. (21) 



We note now that, even in the uniform case Di = 0, ( pT]) has nontrivial periodic solutions 
for Gk"^ — 2X = n'^, with n an integer. Nontrivial solutions for n = are incompatible with the 
periodicity of the nonlocal equation and 6k'^ — 2A = 1 defines the line Fiv, where the constant 
solution bifurcates to solutions of EL that have a sinusoidal perturbation superimposed on the 
constant background. The cases n > 1 define lines which are far from the relevant region, and 
will not be considered. 

There is still another bifurcation line which we know of. We shall see in Section 5 that at 
the line k"^ = j there is a bifurcation from the SC phase. We denote this line by Fni. 

We have verified numerically that in the region between Fni and Fiv there exists a solution 
of EL in addition to (|1^) and ([T5|), which approaches continuously the SC and the constant 
solution at the respective bifurcation lines. This bridging solution is never a minimizer of G. 
The numerical procedure is described in Appendix B. 

The point P at A = /c^ = | is a multiple point that sits at the intersection of four bifurcation 
hues. As a consequence, we expect EL to be highly degenerate at P. To get a better insight 
into the degeneracy, recall that near P one can assume that w is very small. Therefore we 
temporarily drop the cubic term in (^TJ) (cf. below). The modified equation has a three- 
parameter solution at P (for constant D): 



The arbitrary phase 6q is due to the independence of D on 6, the arbitrary normalization A 
is due to the bifurcation at Fi or Fn, and < 7 < 1 is due to the bifurcation at Fni and Fjv 
Since Fni and Fiv cut each other, this family of solutions bridges between SC (7 = 1) and the 
constant solution (7 = 0) at a single point. 

Our range of "interesting temperatures" focuses on the vicinity of P, i.e. where the coher- 
ence length is of the order of the diameter of the ring. 

Let us now outline the case of nonuniform D. If D and w can be expressed by ( p!9D and 
(pOD, then (|2T|) can be solved by Fourier analysis. For conciseness, we take 



wp = A[l 



jcos{9-9o)]. 



(22) 



00 



= ^ Pn COS n9. 



(23) 



n=l 
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Addition of odd terms is straightforward. We make the important observation that if D is 
of type S then Di has the form (p3D, with Pi < 0. li D is not of type S, we still define 
(3i = ^ Jq^ D I cos 9 d9 and take /3i < 0. Unless the first harmonic of Di vanishes, this choice 
of sign is just a choice of the origin of 6. 



The solution of (pif ) is 



71=1 

This expansion provides an approximation for the order parameter in most of the relevant 
region, but is not valid near some of the bifurcation lines. We shall see that, due to the highly 
degenerate character of P, even an infinitesimal deviation from constant D has a profound 
influence on the phase diagram. We have cumulative evidence that, for nonuniform thickness, 
P splits and the phase diagram looks as in Fig. 2. In particular, there are regions where the 
doubly-connected phase resembles more the SC than the constant solution. The following 
sections are devoted to the justification of Fig. 2 to the study of w{9) near the bifurcations. 

We mention that in the phase diagrams (Figs. 1 and 2) we have marked only the two global 
minimizers we know of, singly and doubly connected. We have also met other qualitatively 
distinct solutions of EL, such as the "bridging solution" mentioned above, solutions with 
several oscillations, and solutions with a minimum where D is largest. As far as we have 
tested, these solutions are never global minimizers, and we do not report on them in detail. 

An experimentally interesting quantity is the supercurrent. Far from the bifurcations we 
introduce ( pJ]) and (pOD into (|]) and note that Di and Wi average to zero. We obtain 

/ = khae*{X - e)D/{m*p\) + ©(e'), (25) 

where D, the average of D, is now required. For fixed A, this current attains its maximum 



value at A; = yA/S. Using the definition of A, this maximum is, up to O(e^), 

4. BIFURCATIONS FROM THE NORMAL PHASE (N) 
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When A = and 7^ 0, the global minimizer is obviously y = 0. Increasing A, with fixed 
k, we expect to reach a critical value Ac, corresponding to the critical temperature, at which 
the normal solution will bifurcate to a minimizer yc which does not vanish identically. Here 
the subscript c stands for either I or II. 

In order to perform the bifurcation analysis of ([TTl), we introduce a small parameter rj, and 
expand 

A = Ac + r/A^i..., w = 'qwc + rfwr^2 + ■■■ (27) 
We also write an expansion for A 

+ (28) 

A2 A3 Jo Dwl ' 

with Ac = A(wc) = Jo^{Dwc)~^d6. We thus obtain a sequence of equations for Wc, Wjj2, •••• The 
leading order term Wc obeys 

Dw^w^ + D w,w^ - -iwj' + 2KDwl - -^ij^f = • (29) 

The second order term U7^2 satisfies 

L^(u7^2) = 2X,Dwl - 2Xr,iDwl (30) 

where the linear operator is defined by 

L,{u) = Dwy + {D'w, - Dw'y + {{Dw':) + AKDw,)u - -^-f / ^ . (31) 

We remark that Equation (^) can be written in the form L^{wc) = 0. Notice also that (^9[) 
is invariant under multiplication of Wc by a constant. 

We shall now study the curve Ac = Ac(A;), and the corresponding structure of Wc- Equation 
( P^ is difficult to analyse, as it is nonlocal in a non-standard way. This difficulty can be 
overcome in the case of weakly nonuniform rings. We shall support our findings by solving 
( p9D exactly, with arbitrarily large nonunifomities, for the case where D is piecewise constant. 
The high degeneracy near P implies that two types of expansion are needed. We shall first 
study the bifurcation from the normal state to the doubly connected phase away from P. The 
solution in the vicinity of P is postponed to Section 6. The bifurcation to the singly connected 
phase is a particularly simple case. 

10 



4.1 Bifurcation from Normal to DC 

We are after a nontrivial solution of Equation (^) . We choose the normalization such that 

/■27r 

/ wje = 2n (32) 
Jo 

and consider almost uniform thickness, given by ([T9|) . Using our knowledge for the case of 
constant D, we write Wc = I + ew^i + e^W2 + ... and Xc = k"^ + eXd + e^Xc2 + Expanding 
( p9D and using (^21), gives (to first order in e) 

w'l + Ak^wi = -2Xci - Ak^Di. (33) 



Integrating and applying ( p^ ) and periodicity, we find that Ad = 0. Using ([23|), Fourier 
analysis of (|33D gives 

^ (3n cos ne 

^1 = 1^ -^—r- (34) 

n=l 4/c2 -L 

Except for normalization, this is the limit of (^) for X ^ k"^. 

Expanding to second order in e gives, after some manipulation, 

^ n=l ^ J- 

We see that, no matter how small e is, the approximation scheme breaks down for k"^ = ^. 
This was expected, since at this point Fj intersects other bifurcation lines. 

4.2 Bifurcation from Normal to SC 

As in the previous subsection, Ac is characterized by the fact that (§), or (0) without 
the term involving A, have a nontrivial solution, but the nonlinear term XDy^ or 2XDw^ 
is negligible. It is easier to use the y formulation. Writing yc = Vo + eyi + e^y2 + ... and 
Xc = \ + eXci + e^Ac2 + we obtain a sequence of equations of the form 

Ln{yo) = 0; Lu{yi) = fi, (36) 

where fi is expressed in terms of the functions yj with j < i and 

Liiiyi) = y" + ^Vi . (37) 
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Lii is linear and self adjoint in the space of functions that vanish at and 2tt, with inner 
product {u,v) = Jq^ uvdO. Therefore, fi must fulfill Fredholm's condition for solvability of 
the nonhomogeneous equation: 

rhyode = ^. (38) 
Jo 

This determines Xd- Taking D of type S and recalling the form (^) and the normalization 
521), we obtain 

2/0 = V2sin(^/2), 



3 1°° 



n 



9 ^ Pn ~l~ PnPn—1 

— 1 



(39) 



Note that, since Pi < 0, \c < \- 

4.3 Piecewise Constant Thickness - an Exact Solution 

We note from (0) that the first harmonic diverges for A;^ — s> |. Therefore, wi{0) will be 
a large negative quantity and, for k"^ sufficiently close to |, ewi{0) will not be small and the 
expansion of Subsection 4.1 breaks down. We conjecture that, for arbitrary nonuniformities 
that have a first harmonic component, w{0)/ Jq"" wd9 becomes as A;^ — >• | (along the line 
Fi). To support this conjecture and in order to have a test field for our approximations, we 
consider a special class of nonuniform rings, such that (|29D can be solved exactly. In Section 6 
we shall verify this conjecture for the case of weak arbitrary nonuniformities. 

Let D have the form 

id 9 <tt/2 

D={ (40) 

[ 1 7l/2<e <7T 

with d < 1, and D symmetric about ^ = and about 6 = n. In each interval of the ring D is 
constant, so that differentiating (|29|) and dividing by Wc gives (in each interval) 

wl' + iXcw'^ = 0, (41) 
12 



which has general solution of the form Wc = Ai[l — '~fiCos2y/X^{6 — 6i)]. We are left with 7 
parameters: Ai, 'ji, 6i and Ac- These are determined by normalization, by having the constant 
(27rA;/Ac)^/-D in ( p9l) given by (^) in both intervals, and by continuity of Wc and Dw'^ at 
6 = 71/2 and 6 = 37i/2. Discontinuity of any of these quantities would lead to a divergence of 
Wc or w'^ (with the discontinuity in D regarded as the limit of a fast change). 

Using the symmetry of the problem about and vr, using continuity of Dw'^ and the fact 
that D'^{wcw" + 2XcWc) has the same value at 6' = and 6* = vr, taking Wc{0) < Wc(vr) and 
choosing a simple normalization, we write 

f 1 -7cos2v^e e < 7i/2 

Wc={ ^ ^ ^ (42) 

[ d + d'ycos2y/X'c{'^ - 6) 7r/2 < 6* < vr 
The requirement of continuity of Wc is now 

7cos(7ryA^) = ^-p^. (43) 



Substituting ( ^21) into (|^) and ( ^91) gives, after some trigonometric manipulations. 



7 / /1 + 7 ttVA^, , /I — 7 TTvA^, ,,,, 
TTK = arctan 1 / tan + arctan u / tan . 44 



Equations (B3|) and (33) determine the eigenvalue Ac and Wc{0) = 1—7. The solution that 



describes the line Fj is 

n 2 



Ac 



1 . f2Vd . \ 

— arcsm sm irk 

IT \l + d J 



(45) 



7 = (l-d)/V(l + rf)2-4dsin^7rA; . (46) 

For — i> i, (P6| ) gives 7 — 1, that is, Wc vanishes at 6' = 0, which is the center of the thin 
portion of the ring. This confirms our conjecture. This situation describes the point Pi in 
Fig. 2. 

Writing d = 1 — e, assuming that k is bounded away from |, and expanding in e gives 
Ac — A;^ = — (/i;/47r) tan(7rA;)e^, 7 = e/(2 cos tt/c), in agreement with Subsection 4.1. However, if 
we also write k = ^ — 6 and take the limit e with a fixed ratio 6/e, then 7 = 0(1). More 



precisely, we find 7 = 1/yl + (27r(5/e)^ and Ac — | = — e/(27r7). 
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4.4 The Bifurcating Branch (A > Ac) 



The operator defined in (|3T| ) is self adjoint in the space of periodic functions under the 
inner product 

{u, v) = / ^de. (47) 
Jo 

Therefore, application of Fredholm's condition to (^) gives 

A.i = \c r Dwlde/ Dw^O- (48) 
JO Jo 



This determines the value of 7] in (^Tf). It is easy to show that (|48| ) holds for the SC solution, 
too. 



Substituting now (^) into (|T2|) we find that to leading order the free energy is 

AttXJo^ Dw^de 

5. BIFURCATIONS FROM THE SINGLY CONNECTED PHASE 



(49) 



The present section deals with an essential property of our model, and the main results 
are valid for arbitrary shape of D{6), i.e., the nonuniformity in the thickness is not necessarily 
small. For simplicity, we shall assume again that D is of type S. 

In the SC solution w{0) = and near a bifurcation from it w{0) is small. The existence and 
nature of the bifurcations from the SC solution stem from the divergence of A as w{0) 0. 

Let us first assume that D{6) is smooth in the region near 6 = 0. Let us write for shorthand 

fiO) = D{9)w{e), ~f{9) = /(O) + ln0)9', (50) 

where we have used the symmetry of the problem to include only even powers of 6, and let us 
separate A into a divergent and a regular part: 



with 



2 J. f Jo ff 
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Ar remains finite for w{0) 0. Performing the integration in 



TT 



A = + Ar. (53) 



We set now 6* = in (IT^), note that w'{0) = 0, introduce ( |53D and (|50D, and expand to 
order w (0)^/^. This gives 



^2w{0)w"{0)D{0)Ab, = 7r(l - {2k)-^). (54) 

The bifurcation hne is obtained by taking w{0) 0. This gives k = ^, independently of A. 

If Ar > (resp. Ar < 0), then (^) can be satisfied with w{0) > only for k"^ > j (resp. 
fc^ < |). We can therefore expect that for Ar < there is a bifurcation from the SC to the DC 
phase and, for Ar > 0, to the bridging solution in the irrelevant region. Numerical solution 
of {\n\ } for fc^ ^ i and various values of A confirms this scenario. 

If the thickness is almost uniform, w{6) in the SC phase can be approximated by means of 
( P^D and (ppp. Once w is known, Ar can be evaluated. Keeping one correction term in e and 
in 7], we obtain 

Recalling that Pi < 0, we see that Ar < for A sufficiently close to Ac, but eventualy changes 
sign at a point which we may identify as P2 in Fig. 2. Taking Ac from (^), we find that P2 is 
located at A = ;| — e/3i/2. 

If the thickness is uniform, then Ar = at P and is positive for A > Ac, in agreement with 
Fig. 1. 

Let us finally discuss the case where D'{6) is discontinuous at = 0, as in the example of 
Ref. 1^. Due to the continuity of Dw', the discontinuity of f'{0) will be of order w{0). The 
same will be true for the coefficient of 6^ in the expansion of /. Therefore, (^) can still be 
used and (|52|) does not diverge as w{0) 0. 

6. THE REGION NEAR THE CRITICAL POINTS 
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We shall again restrict ourselves to the case where D is almost uniform (cf. (|T9D), define 
5 = \ — k"^ and Ai = (A — |)/e, and expand around the point P: 

w = Tjiwp + ewi + ...). (56) 

The expansion is in powers of e only, but we explore the region where Ai, 5/e and rj/e are 
finite and the terms in the expansion depend on these parameters. The leading order term wp 
is given by the expression (^), where without loss of generality we set A = 1 and the other 
parameters are still undetermined. At this stage it will be instructive to consider D which is 
not necessarily symmetric about and vr. 

Substituting (|56|) into ( pTD leads to the equation Lp{wi) = fi, with 



Lp{u) = wpu —WpU +u — / -^d6, (57) 

ZTT JO Wp 



u 



= _ADi - wpWpD, - 2\iwl + ^wl + —— / ^de - 2A-, (58) 

2e ^ ZTT Jo Wp e 



' / — t 

2n Jo Wp 

and A = 1 — 7^. 

As in the situations considered above, Lp is linear and self adjoint under the inner product 
{u,v) = {uv / wy)d9 . Also, it is easy to check that Lp{wp) = 0. The unusual feature is that 
there are other solutions as well, that is, Ker(Lp) is not one dimensional. In fact, Ker (Lp) 
is three dimensional. It is spanned by the three functions {l,wp,Wp}, corresponding to the 
three parameters m. wp. 

Therefore we obtain a set of three solvability conditions 



r4de=r^de=r^de = o. (59) 

JO Wp Jo Wp Jo Wp 



The last condition implies 

Di sm{9 - 9o)d9 = 0. (60) 

JO 

This equation determines the position 6 = 6q a.t which wp has a minimum and may eventually 
vanish. Note that 6q is determined by the geometry only, and does not vary with the tem- 
perature or the magnetic field. After having obtained (160|) , we can redefine 6 and set 6q = 
without loss of generality. 
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The other two conditions in (^) reduce to 



and 



Ai + A-i/2^ = ^ (61) 



4Ai + 4AV2^ - - (1 + = 0. (62) 

e 2 e 



These two equations determine 7 and t]. The following formula that can be derived from 
(^l])-(^) will be useful later on: 

= (63) 

e 3 37 

6.1 The Bifurcation Line Ti Near the Point Pi 

This is obtained by setting r/ = in ( |5^ and (pT]). The lowest value of A for which both 
equations are fulfilled is obtained for 

Ai = A,i = -ypj + {A6/ey, (64) 

and 

/3i /3i .... 

7 = , = TT— • (65) 

^(3! + (45/e)2 4A,i 



The limit considered at the end of Subsection 4.3 is a special case of ([65|), with /3i = — 2/7r. 
It is also interesting to note that (|6^ ) can be extrapolated to the case that 5 is not small and 



it still reproduces to order e the result obtained for Ac in Subsection 4.1. 

We note that no matter how small e is, i.e. even for arbitrarily weak nonuniformities, 
7 — s> 1 as 5 — » 0. Hence we deal here with a pattern selection mechanism, where the SC 
pattern is selected among all the infinitely many solutions (^2]). 

6.2 Evaluation of wi 



For conciseness, we shall again consider D of the form (|23D. Fourier analysis of the equation 
Lpiwi) = fi and use of (|6^) yield 

00 

Wi = ^ a„ cos n9 



n=0 
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with 



do = 1 

^ 3 3 



and 



^(/5„_i-/5„+i) + 3f^ (66) 



for n > 3. The coefficients ao and ai remain undetermined to the present order in the expansion 
(^6]); this is the same situation we met in Eq. (|2^) , where A and 7 remained undetermined. 

6.3 The Region Near the Point P2 

We know from Section 5 that near (fc^ = = i ~^/^i/2) there is a special point. The goal 
of this subsection is to obtain the behavior of w(0) ~ riwp{Q) near this point. For this purpose 
we eliminate rj from (^) and (0) and substitute 7 = 1 — wp(0). This gives an equation for 
W7p(0): 



e 3 1 — wp[y)) 2 
Near P2, wp{0) « 1, hence it is negligible compared to 1 or 2. However, in the relevant 
regime \iWp{Q) is of the same order as Ai + /3i/2. Therefore, (|67|) becomes a third degree 
equation for wp(0)^/^. The solution of this equation provides the missing information about 
'w{6) to present order in the region near P2. 

We see from (^) that the relevant region \ — k"^ = 5 > Q corresponds to lyp(O) > 1 + 



/9i/(2Ai). In order to characterize the solutions of (|67|), we have drawn in Fig. 3 a contour 
plot of (5(wp(0), Ai). At the axis wp(0) = 0, 5 = 0. For Ai < -/9i/2, all the factors in (|3) 
are positive, 5 increases monotonically with wp(0), and there is only one wp(0) for a given 
5. However, for Ai > — /3i/2, there is initially a valley into the irrelevant region, that is, 5 
is negative - a situation which cannot correspond to the minimal energy. The level 5 = is 
reached again for wp(0) = 1 + /3i/(2Ai), and this is found to be the global minimizer. This 
situation corresponds to the ray to the right of P2 in Fig. 2. The line Fiv is characterized by 



a<5 



0; this happens at 5/e ^ -(l/(3v^))(2(2Ai + Pi)/3f/^ 



dwp(0) 

6.4 Higher Order Treatment 
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The first order expansion in ( p6D leaves the O*'^ and l'^* harmonic of w{6) undetermined 
for the doubly connected phase. This may not be a grave shortcoming, unless they diverge 
and their contribution is actually of 0(1). We might especially suspect that the behavior of 
the entire w{0) near P2 is qualitatively different from that of wp{0). We therefore proceed to 
evaluate these harmonics near P2. 

In order to keep the analysis simple, we now restrict ourselves to D{6) of the form D = 
l + e/?i cos 6' with {f3i < 0). Cases of nonuniformities that are neither smooth nor small will be 
reported elsewhere. Using (|56D, ( ]63| ) and ( ]66| ) sufficiently close to /c^ = 1/4, Ai = —Pi/ 2 and 
recalling that 7 = 1 at P2, we can write 

w = -2pie{l+pe-cos9){l + qe- Piecose) + 0{e^), (68) 

with p = ao + ai + Pi/2 and q = —ai — Pi. Our strategy is now to use (|68|) without the 
correction term as a variational form, and to find the parameters p and q that minimize G in 
([Tol). Since w{0) is proportional to p, we are especially interested in the behavior of p{\i,6) 
at this minimum. 

Eq. (|T^) contains two integrals that require careful handling. One of them is / D{w'Y /wdO, 
which contains a vanishing denominator for pe 0. In this case we perform the integral with- 
out taking pe as a small quantity, and expansion on it is performed only after the integration. 

The second delicate integral is A. Using the periodicity and the behavior of cos6' in the 
complex plane, we can write A = J{wD)^^d6, where the integral is around the rectangle with 
vertices at — vr/2, 3tt/2, 3tt/2 + ooi and — 7r/2 + ooi. Since wD is the product of three simple 
factors, it is easy to find the poles and residues of the integrand, and A can be evaluated 
exactly. 

We write now A;^ = ^ — 5, with <C< 00, A = 1/4 — ePi/2 + e^A, and expand G up 

to 0{e^). We obtain a polynomial in and q. Minimizing this expression gives 

8A + 2/5i(2j9 + 3g) -9/^1 = (69) 

and 

^ = ^p/(-yP-A + — ). (70) 
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One immediately recognizes that ([70|) describes the same behavior as (0), but the point P2 
is shifted to 

Ap. = ^-^6A + ^(eA)'. (71) 
7. DISCUSSION 

We have analytically obtained the phase diagram in the entire temperature-magnetic flux 
plane, as schematically drawn in Fig. 2. Some of our results are valid for small deviations 
from uniform thickness, some results are valid for arbitrary shape of D{9) and some are valid 
for piecewise constant thickness; nevertheless, all the results are consistent among themselves 
and with Fig. 2. In addition, we have numerically tested the predicted dependence of the 
minimum value w(0) of the order parameter as the physical parameters A and k are varied. 
An integrative look of all the cases considered makes it reasonable to expect that the results 
found here will remain qualitatively correct for almost any form of the thickness D{9). 

The results obtained here enable us to check the numerical flndings obtained in Ref. |^ 
for k"^ = ^ and for a particular form of D{6). In we compared the free energies of the SC 
and the DC phases. Upper bounds for these free energies were obtained by using sequences of 
Fade approximants of increasing length as variational forms. G in equation (|^) was minimized 
in the space of the coefficients of these approximants. In 0] we did not suspect that the DC 
phase goes continuously into the SC phase at the line Fm. Since at this line both phases 
become the same, our previous approach was clearly inappropriate. Indeed, there is a region 
in which the singly-connected phase is thermodynamically favorable. We also confirm here 
that along the line k"^ = ^ there are two phase transitions as the temperature is varied below 
Tc. However, the transition at F2 is not discontinuous (first order), as concluded from our 
early numerical results, but rather a continuous (second order) transition after which w{0) 
increases at the rate dw{0)/dX = 4. 

In the relevant region (A;^ < j) the stability domain of the singly-connected phase is 
restricted to a line segment, and has measure zero. This property might have applications for 
calibration purposes. However, the transition to the doubly-connected phase is continuous, 
and near this segment w{0) will be small. 
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A simple experimental verification of our predictions could be performed by measuring the 
supercurrent by means of a magnetometer, as in Ref . . Fig. 4 shows numerical evaluations of 
the supercurrent, given by as functions of the magnetic flux for three different temperatures 
close to P2. 

Another experimental verification of our predictions could be obtained by detecting the 
current fluctuations 0. In the doubly-connected phase, for A > Apj, the supercurrent does 
not vanish when the flux is a half-integer number of quanta (fc^ — \)- As a consequence, 
this situation is degenerate, since the supercurrent can flow in either sense. Near P2, the 
energy barrier between both minima of G is small, and the supercurrent fluctuates. However, 
for A < Ap2 the supercurrent vanishes, and fluctuations should disappear, while the sample 
remains superconducting. 

A third way to detect whether G has one or two minima is by the absence or presence of 
hysteresis as the magnetic field is varied: for A > Apj, if the flow $ is slowly raised and crosses 
a half integer of quanta $0; then according to (^) and (|^, k ought to jump from | to — | 
in order to minimize the free energy. Following and (0), this implies a jump in the order 
parameter. Since the order parameter stays in a local minimum of the free energy and the 
sample 'does not know' that there is a still lower minimum, the jump in the order parameter 
and the inversion of I may occur only after k is larger than |. If the change in $ is now 
reversed, the order parameter will jump back at a lower value of $, after k is smaller than |. 
We argue that this hysteresis phenomenon should not occur for A < Apj. Although at = ±| 
there is still a jump in the phase 0, equation (||) suggests that this jump concentrates where 
u vanishes, and is thus not physically observable. 

Our analysis for almost uniform thickness shows that the important part of the nonuni- 
formity is its first harmonic. This might explain why the effects predicted here were not 



observed in [jTO| for samples with a weak link (a small region where D(6) « 1), which might 
be thought to be extremely nonuniform. When the deviation from uniformity is not small, 
Eq. (p9D suggests that coupling between different harmonics becomes important. 

From a mathematical point of view, the the Little-Parks problem poses a nonstandard 



21 



minimization problem: 

Min G[y]k]= r iy'fDde + i2T:kfK-^ (72) 
io 

subject to the constraint Jq^ y^DdO = 1. It might appear that this is a very specific problem, 
but it actually stems from very general principles: the Ginzburg-Landau equation (which is 
nothing but a slowly varying small order parameter expansion), conservation of current and 
single- valuedness of the order parameter. In general, for a multiply connected sample, we shall 
have a problem with a larger parameter-space. Instead of one k we shall have a ki for each 
hole in the sample that can enclose a flux and we shall have a Aj for the supercurrent around 
that hole. 
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Appendix A. Influence of the Expelled Flux 

We deal now with the contribution of the term {H — HextY in dH), which has so far been 
neglected. Let us first consider the case where the sample is a long circular cylinder. Then 
it is well known that H — H^xt = 4:TTl/{hc), where h » R is the height of the cylinder. 
This additional magnetic field exists inside the cylinder, i.e., in a volume irR^h. Therefore, its 
contribution to G is 

Gh = 4kVK-'^h-^\A-y Dd9, (73) 

Jo 

where we have used @, R = 1, and k = Xl/^ is the Ginzburg-Landau parameter. Note that 
Gh is proportional to the thickness D/h, and we are no longer free to choose the normalization 
ofD. 

If the sample is not a long cylinder, then Gh has to be divided by the self-inductance of the 
cylinder and multiplied by the actual self- inductance of the sample. For simplicity, let us stick 
to (^). The infiuence of Gh in the EL equation is that the nonlocal term in ( pT]) is multiplied 
by the factor (1 + 27iK^^h^^XA^^). Therefore, a sufficient condition for this correction to be 
uninportant is K.~'^\D/h « R. 

We are mainly interested in the inffuence of Gh when A diverges. Following the steps of 
Section 5, we reach again (|5^ , but A/j has to be replaced by Kr — ttA/ {hn?). Therefore, the 
bifurcation from SC is still along the line k'^ = but the point P2 moves to higher values of 
A. 

Appendix B. Numerical Solution of Equation ( [Tl] ) 

We used the built-in MATHEMATICA functions. The procedure is as follows. 

1. Fix ti;(0), A and C = 2{2'nk/Kf. 

2. With these values fixed, we can use NDSolve to evaluate w{9). This function requires 
the initial values of w. w{Q) is given and w'{Q) = 0. One has to be careful not to pick 
values of C which are too large, since they lead to w(7r) >> 1. 
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3. The true value of C is now obtained by requiring periodicity. If w is expected to be of type 
S, we can use the function FindRoot to require, for instance, that w(3) — w(27r — 3) = 0. 
(This difference is now regarded clS cl function of C.) 

4. After the true C and the corresponding w{9) are known, we obtain A by numeric inte- 
gration. We then evaluate /c^ = CA^/(87r^). In this way we obtain k"^ as a function of A 
and iu(0). 

Figure Bl shows a sequence of solutions of equation (|n|) with varying values of w{Q). 
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Figure Captions 



Fig. 1: Phase diagram for a ring of uniform thickness. A is proportional to Tc — T 
and k measures the deviation from an integer number of magnetic flux quanta. The dashed 
lines denote bifurcations between solutions of the Euler-Lagrange equation which are not 
global minimizers. The line with uneven dashing denotes first order (discontinuous) transition 
between the singly and the doubly connected phases. 

Fig. 2: Phase diagram when the thickness is not uniform. P splits into two critical points, 
Pi and P2, and the portion of the the bifurcation line Pni between these points is physically 
relevant. 

Fig. 3: Contour plot of 5 = | — A;^. 6 increases toward the lower right. At the upper left 
there is a "valley" , where 6 < 0. The line Piv in Fig. 2 "goes down" along this "valley" . 

Fig. 4: Supercurrent as a function of k for several values of the parameter A = {R/C,)"^. 
For A < Ap2, the supercurrent vanishes at = |. For A = Apj, the supercurrent rises with 
infinite slope. For larger A, /(|) 7^ 0. /max is given by Eq. (p6l). 

Fig. Bl: w{6) and k as functions of w{0). The value obtained for k is shown next to the 
corresponding curve. For the present illustration we took uniform D and A = 1. This is the 
"bridging solution" , discussed in Section 3. The extreme situations were taken from ( [T3| ) and 
(pISl). If D is weakly nonuniform and k ^ 0, then the DC phase looks qualitatively as this 
bridging solution. 
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